Abstract 

Piloting difficulties associated with conducting 
maneuvers in hypersonic flight are caused in part 
by the nonintuitive nature of the aircraft response 
and the stringent constraints anticipated on allow- 
able angle-of-attack and dynamic pressure variations. 
This report documents an approach that provides 
precise, coordinated maneuver control during excur- 
sions from a hypersonic cruise flight path and ob- 
serves the necessary flight condition constraints. The 
approach is to achieve specified guidance commands 
by resolving altitude and cross-range errors into a 
load factor and bank angle command by using a co- 
ordinate transformation that acts as an interface be- 
tween outer- and inner-loop flight controls. This in- 
terface, referred to as a “resolver,” applies constraints 
on angle-of-attack and dynamic pressure perturba- 
tions while prioritizing altitude regulation over cross 
range. An unpiloted test simulation, in which the re- 
solver was used to drive inner-loop flight controls, 
produced time histories of responses to guidance 
commands and atmospheric disturbances at Mach 
numbers of 6, 10, 15, and 20. These time histories are 
used to illustrate the manner in which the overall con- 
trol system incorporates flight condition constraints 
and accounts for high-speed flight effects. Angle-of- 
attack and throttle perturbation constraints, com- 
bined with high-speed flight effects and the desire 
to maintain constant dynamic pressure, significantly 
impact the maneuver envelope for a hypersonic 
vehicle. Turn-rate, climb-rate, and descent-rate lim- 
its can be determined from these constraints. Den- 
sity variation with altitude strongly influences climb- 
and descent-rate limits and throttle modulation if 
dynamic pressure is to be maintained during verti- 
cal transitions between cruise flight conditions. 

Introduction 

Many recently proposed hypersonic vehicle con- 
cepts present significant problems in guidance, flight 
control, and flying qualities that must be addressed 
if such designs are to be made viable. These air- 
breathing single-stage-to-orbit (SSTO) vehicle de- 
signs traverse a broader range of flight regimes than 
aircraft that have flown in the past, and must em- 
phasize performance during all phases of flight from 
takeoff to orbit to achieve their mission. Stringent 
constraints on angle-of-attack, sideslip, and dynamic 
pressure variations must be observed, since hyper- 
sonic propulsion system performance is strongly de- 
pendent on flight condition and may suffer dramat- 
ically if excessive variation is experienced. Also, an 
unprecedented level of coupling exists between flight 
dynamic, propulsion, and structural modes, and the 


degree of this coupling may vary significantly over 
such a large design envelope. The anticipated sensi- 
tivity of the propulsion system to variations in flight 
conditions makes precise regulation of angle of at- 
tack, angle of sideslip, dynamic pressure, and flight 
path imperative, while the high degree of coupling 
makes such tight control difficult to achieve. The 
area of maneuver control for hypersonic vehicles re- 
quires further investigation to reveal the various pi- 
loting issues and to determine the level of automation 
necessary to achieve the desired performance. 

Previous work by Berry (ref. 1) involved a piloted 
simulation of a hypersonic configuration in which nu- 
merous issues in handling qualities relevant to ma- 
neuvering in hypersonic flight were examined. Recent 
research by Lallman and Raney (refs. 2 to 6) has been 
focused on regulation of trajectory parameters dur- 
ing hypersonic maneuvers at cruise flight conditions, 
but maneuvers were restricted to the vertical plane. 
This report is focused on the problem of coordinated 
control of maneuvers in both the vertical and lat- 
eral planes at hypersonic cruise flight conditions. An 
approach is presented that provides maneuver coor- 
dination through the resolution of altitude and cross- 
range errors into a combination of normal load factor 
and bank angle command. The goal of this research is 
to define an automatic control design option for ex- 
ecuting coordinated maneuvers in hypersonic flight 
while regulating key parameters, such as angle of at- 
tack, angle of sideslip, and dynamic pressure. Addi- 
tionally, this work is intended to be an investigation 
of the manner in which propulsion and flight condi- 
tion constraints impact the maneuvering capabilities 
of hypersonic aerospacecraft. 

A control integration concept is described in this 
report for powered hypersonic cruise-turn maneu- 
vers. First, several introductory remarks on critical 
issues relevant to the control of air-breathing vehi- 
cles in hypersonic flight are made. An overview of 
the basic components of a conceptual control system 
for an air-breathing hypersonic vehicle is then pre- 
sented. The overall control system is divided into 
three main subsections consisting of (1) inner-loop 
controls to provide stability augmentation, (2) outer- 
loop controls to track guidance commands and reject 
disturbances, and (3) an interface between the inner 
and outer loops that resolves vertical and lateral ac- 
celeration commands into a lift vector command that 
is specified by a normal load factor and bank angle 
combination. Each of these subsections is described 
in order following the control system overview. 

The interface between inner and outer loops, re- 
ferred to throughout this report as a “resolver,” is the 
main subject of this research. The resolver concept 



is introduced following the discussion of the outer- 
loop controls. Logic associated with the resolver to 
prioritize regulated trajectory parameters and incor- 
porate flight condition constraints is then developed. 
Following this description of the resolver, a throttle 
control law, which regulates dynamic pressure during 
hypersonic maneuvers, is presented. The manner in 
which constraints on angle-of-attack and thrust per- 
turbations impact the maneuver envelope of a hyper- 
sonic vehicle is also discussed. A control design case 
is presented for an example vehicle, and a test simula- 
tion is described which was used to produce time his- 
tories of responses to guidance commands and atmo- 
spheric disturbances. These responses illustrate the 
manner in which the combination of inner-loop, re- 
solver, and outer-loop control systems operate. The 
simulated responses are discussed, and several con- 


elusions are presented that have significant implica- 
tions for maneuver control during hypersonic flight. 

Symbols 

Ay 

lateral acceleration in body-axis 
coordinate system, ft/sec 2 
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wing span, ft 

C D 

drag coefficient 
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top 

Cl 

lift coefficient 

Cl, a 

_ dCLp 
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C p rate of density variation with 

altitude, slugs/ft 3 /ft 

D drag force, lb 


g acceleration due to gravity, 

ft/sec 2 

h altitude, ft 

I x rolling moment of inertia, 

slugs-ft 2 

I xz product of inertia, slugs-ft 2 

I z yawing moment of inertia, 

slugs-ft 2 

Kap sideslip feedback gain to aileron 

Kag r rudder deflection cross-feed gain 

to aileron 

Kap bank angle feedback gain to 

aileron 


Ka j 

<p 

roll-rate feedback gain to aileron 

Krp 

sideslip feedback gain to rudder 

Kr r 

yaw-rate feedback gain to rudder 

I<h 

feedback gain for proportional 
altitude error 

K hd 

feedback gain for derivative 
altitude error 

K h i 

feedback gain for integral altitude 
error 

K- q 

feedback gain for dynamic 
pressure error 

Ky 

feedback gain for proportional 
cross-range error 

K yd 

feedback gain for derivative cross- 
range error 

Kyi 

feedback gain for integral cross- 
range error 

L,M,N 

rolling, pitching, and yawing 
moments, ft-lb 

Mach 

ratio of vehicle velocity (relative 
to air mass) to speed of sound 

m 

mass, slugs 

n 

normal load factor with respect 
to vehicle flight path, g units 
(1 g = 32.17 ft/sec 2 ) 

n t 

lateral component of normal load 
factor with respect to an Earth- 
fixed axis system, g units 

n v 

vertical component of normal 
load factor with respect to an 
Earth-fixed axis system, g units 

P,Q,r 

roll, pitch, and yaw rates in 
body-axis coordinate system, 
deg/sec 

Q 

dynamic pressure, lb/ft 2 

R 

altitude measured to center of 
Earth, ft 

S 

wing area, ft 2 

s 

Laplace transform parameter 

T 

thrust force, lb 

U, V, w 

body- axis velocities, ft /sec 

V 

inertial velocity, ft/sec 
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X, Y, Z 
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V> 


Ll> 


weight, lb 

force components of body 
axes, lb 

vehicle body axes 

cross range, ft 

angle of attack, deg 

sideslip angle, deg 

perturbation value or deviation 
from nominal 

aileron deflection, positive right 
trailing edge down, deg 

rudder deflection, deg 

damping ratio 

pitch angle, deg 

distance from center of gravity 
to lateral accelerometer, positive 
forward, ft 

air density, slugs/ft 3 

real root in complex plane 

complementary filter time 
constant 

roll angle, deg 

heading angle, deg 

natural frequency, rad/sec 


Abbreviations: 

APAS aerodynamic preliminary analysis 

system 

B-T-D bank-to-dive 

p-i-d proportional-integral-derivative 

POST program to optimize simulated 

trajectories 


a parameter used in aileron control 

law 

ace accelerometer signal 

cmd commanded value 

err error signal 

h parameter used in altitude 

regulator 

max maximum allowable value 

min minimum allowable value 

prev value from previous time step 

r parameter used in rudder control 

law 

t parameter used in throttle 

control law 

y parameter used in cross-range 

regulator 

0 nominal or trim value 


Dimensional stability and control derivatives: 
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The coordinate-system notation and reference frames 
used in this report are shown in figure f. 


Specialized terms: 


Superscripts: 

first derivative with respect to 
time 

second derivative with respect to 
time 


altitude priority: the practice of assigning priority to 
the elimination of altitude errors over the elimi- 
nation of cross-range errors 

ballistic origin: origin of resolver coordinate system 
that designates the condition at which normal 
load factor is zero 


~ estimated value 

~ complementary filtered value 

Subscripts: 


bank-to-dive: the practice of rolling the lift vector 
off of vertical to achieve the desired descent rate 
while remaining above some minimum angle of 
attack 
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centrifugal relief: reduction in normal load factor 

required to maintain constant altitude at hyper- 
sonic speeds that results from the centripetal ac- 
celeration of the vehicle as it circles the Earth 

complementary filter: a low-pass filter that is aug- 
mented with derivative feedback in such a way 
as to eliminate its effect on closed-loop stabil- 
ity while attenuating responses to atmospheric 
disturbances 

resolver: a rectangular-to-polar coordinate transfor- 
mation that converts vertical and lateral accelera- 
tion commands into a load factor and bank angle 
command vector while applying limits and logic 
to incorporate features that limit altitude prior- 
ity, “bank-to-dive,” and load factor 

Issues Relevant to Hypersonic Flight 
Control 

Several factors complicate maneuver control for 
air-breathing vehicles at hypersonic speeds. For ex- 
ample, strict regulation of dynamic pressure may be 
required to meet the constraints of a given ascent 
trajectory. This practice is particularly critical in 
the case of air-breathing SSTO vehicles, which may 
follow a trajectory profile determined by a combina- 
tion of engine performance and heating criteria. At 
hypersonic speeds, dynamic pressure is strongly in- 
fluenced by the variation of air density with altitude. 
Compensation for this effect during vertical maneu- 
vers will almost certainly require some form of au- 
tomatic feedback regulator, since it would otherwise 
represent a formidable increase in manual work load. 

Also, lateral maneuvering in hypersonic flight is 
complicated by a high-speed flight effect that man- 
ifests itself as a reduction in load factor required to 
maintain constant altitude. A vehicle in hypersonic 
flight may approach orbital velocities while still in the 
atmosphere. The centripetal acceleration of the ve- 
hicle as it circles the Earth then becomes significant 
enough to cause a noticeable reduction in the load 
factor required to maintain level flight. This effect, 
referred to as centrifugal relief, can have an apprecia- 
ble impact on the bank angle required to complete a 
given lateral maneuver. For instance, a conventional 
aircraft in level subsonic flight that is executing a co- 
ordinated turn with a bank angle of 60° experiences 
a normal acceleration of 2 g units. However, a vehicle 
in hypersonic flight with a velocity of f8 500 ft/sec 
(about Mach 18 at 95 000 ft) is required to bank ap- 
proximately 75° to produce a 2 g turn (ref. 7). In 
straight and level flight about a spherical Earth, the 
normal load factor n required to maintain constant 


altitude is 


where V is the inertial velocity of the vehicle, R is 
its distance from the center of the Earth, and g is 
acceleration due to gravity. In subsonic flight, the 
term V 2 / Rg is negligible, but at hypersonic speeds 
the term becomes significant. The influence of this 
term becomes more pronounced as the velocity of the 
vehicle approaches orbital speeds. 

Another factor in hypersonic flight is the ex- 
tremely slow response of trajectory variables such as 
heading or flight-path angle. Relatively long time 
periods are required to redirect the velocity vector, 
since the momentum of the vehicle is so great. Some 
form of predictive flight director may be necessary 
to guide the pilot through sustained and somewhat 
nonintuitive maneuvers. In addition to the compli- 
cations that are inherent in maneuvering in the hy- 
personic regime, the current military specification 
manuals (Mil Specs) cannot provide adequate de- 
sign criteria for the synthesis of a fly-by-wire system 
for hypersonic air-breathing vehicles. Shortcomings 
have been identified in the areas of roll control for 
hypersonic turns, response to atmospheric perturba- 
tions, cockpit displays for hypersonic maneuvering, 
and many other critical flight-control criteria (refs. 7 
to 9). Also, the exact role of the pilot during the as- 
cent portion of a transatmospheric flight is uncertain. 
Hypersonic propulsion modules are likely to be sensi- 
tive to angle-of-attack and angle-of-sideslip perturba- 
tions, which may disturb inlet conditions and result 
in poor engine performance (ref. 10). Such transients 
in inlet conditions could also result from responses 
to abrupt guidance commands. Extremely low tol- 
erance to perturbations from the nominal flight con- 
dition may require that some tasks be automated if 
the desired mission is to be achieved. 

The control design presented in this paper ad- 
dresses several issues that are of particular signifi- 
cance to maneuvering in hypersonic flight. The de- 
sign concept presented could be used to control the 
vehicle directly or to drive a cockpit-display flight 
director based on inputs from an onboard guidance 
algorithm. The concept acts as an interface be- 
tween the onboard guidance system and the inner- 
loop flight controls. 

An overview of the basic components of a con- 
ceptual control system for an air-breathing hyper- 
sonic vehicle follows, after which each subsection of 
the control system is described. Following the con- 
trol system description, a design case is demonstrated 
for an example vehicle, and a test simulation is pre- 
sented. Simulated responses to guidance commands 
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and atmospheric disturbances are discussed and are 
followed by several concluding remarks relevant to 
the control of air-breathing vehicles in hypersonic 
flight. 

Control System Overview 

The basic elements of the hypersonic control sys- 
tem concept are presented in figure 2. The over- 
all control system is divided into three main sub- 
sections: (f) inner-loop controls to provide stability 
augmentation, (2) outer-loop controls to track guid- 
ance commands and reject disturbances, and (3) an 
interface between the inner and outer loops, which 
resolves vertical and lateral acceleration commands 
into a normal load factor and bank angle command 
vector. A summary of the overall system is presented 
in this section, and each element is described in detail 
in subsequent sections. 

An onboard guidance system is envisioned that 
will issue commands in terms of altitude and cross 
range to achieve a given task, such as executing a 
heading change or eliminating vertical and lateral 
offsets from a desired trajectory profile. These com- 
mands are differenced from actual feedback values, 
as shown in figure 2, to produce error signals. Al- 
titude and cross-range control laws that translate 
errors from the desired trajectory parameters into 
vertical and lateral acceleration commands are then 
applied. This feedback scheme makes up the outer- 
loop control system, which consists of two sepa- 
rate regulator loops to track altitude and cross-range 
guidance commands. The specific architectures of 
the outer-loop altitude and cross-range control laws 
are described subsequently in this report. 

Vertical and lateral acceleration commands from 
the outer loops are converted into load factor compo- 
nents that take into account the reduction in vertical 
load factor that is caused by centrifugal relief. These 
load factor commands are then fed into a rectangular- 
to-polar coordinate transformation, which converts 
the vertical and lateral load factor components into 
a normal load factor and bank angle command vector 
(fig. 2). Logic and limiting strategies are associated 
with the coordinate transformation to incorporate 
various flight condition constraints and to prioritize 
the regulation of trajectory parameters. The coor- 
dinate transformation and its associated limits and 
logic are referred to as a resolver throughout this re- 
port. This resolver acts as an interface between the 
outer loops, which track guidance commands, and 
inner loops, which provide stability augmentation. 

In this study, the load factor and bank angle 
command vectors produced by the resolver were fed 


directly into an inner-loop flight controller. (In a 
real-time investigation, the outputs from the resolver 
might have been used to drive a cockpit display to 
aid the pilot in conducting the desired maneuver.) A 
simple inner-loop controller is presented in this re- 
port to allow simulation of the completed system, al- 
though more sophisticated inner-loop designs can be 
used with the resolver strategy. The inner-loop con- 
troller commands rudder deflections, elevon deflec- 
tions, and thrust variations to track the load factor 
and bank angle command vectors from the resolver 
while maintaining constant dynamic pressure during 
a maneuver. The inner-loop control design is param- 
eterized in terms of specified frequency and damping 
for the desired response characteristic and requires 
values for the aerodynamic stability and control co- 
efficients of the vehicle. 

Inner-Loop Flight Controls 

A simple inner-loop controller was developed to 
allow simulation of the overall control concept. Only 
lateral flight-control laws were designed, and per- 
fect pitch control was assumed for the test simula- 
tion. Short-period dynamics and propulsive thrust 
and moment variations with angle of attack were 
not considered in this inner-loop design. Inclusion of 
these effects would not invalidate the resolver design 
concept presented in this report, but it would neces- 
sitate a more elaborate inner-loop controls design. 
The inner-loop controller acts to stabilize the lateral 
flight dynamic modes, which provides acceptable fre- 
quency and damping characteristics while following 
commands issued by the resolver interface with the 
outer loop. It is crucial that the angle of sideslip be 
kept to nearly zero during any maneuvers, so that 
engine inlet tolerances are not violated and so that 
excessive drag is not produced (i.e. , the maneuvers 
should be coordinated). 

The approach used in this study has been to 
develop two classical feedback loops to drive the 
ailerons and rudder. The aileron control loop tracks 
bank angle commands that are issued by the re- 
solver, while the rudder control loop maintains zero 
sideslip throughout the maneuver. The effect of rud- 
der deflection on rolling motion is compensated by a 
cross-feed to the aileron control law. A simple pole- 
placement technique is used in each of the two control 
loops to determine the required feedback gains. De- 
sired inner-loop dynamics are specified in terms of 
frequency and damping, which allows the designer to 
place the closed-loop poles at a location that might 
be desirable from a handling-qualities perspective. 
Derivations of the rudder and aileron control laws 
are presented in appendix A. 
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The feedbacks used in the rudder control law 
are yaw rate and sideslip angle, while those used 
in the aileron control law are roll rate, bank angle, 
and sideslip angle. Angle of sideslip is an impor- 
tant measurement that appears in both these con- 
trol laws. However, the reliability of air-data mea- 
surements during some phases of hypersonic flight 
is uncertain. In some cases, it may be necessary 
to augment traditional air-data sensors with an es- 
timate of air-data measurements based on inertial 
instruments. For demonstration purposes, the inner- 
loop control system designed in this report utilizes an 
angle-of-sideslip feedback that is constructed from a 
lateral-accelerometer reading rather than an air-data 
sensor array. The lateral accelerometer is placed at 
the instantaneous center of rotation in response to a 
rudder deflection; this placement eliminates the need 
to correct the accelerometer reading for rudder ef- 
fects. A detailed description of the angle-of-sideslip 
estimator is included with the formulation of inner- 
loop controls in appendix A. The use of air-data esti- 
mates based on the designer’s knowledge of the aero- 
dynamic model makes the control system sensitive 
to the fidelity of the aerodynamic predictions. Such 
air-data estimates may be included in a weighting 
scheme in which measurements from various sensors 
are combined to derive a best estimate of the quan- 
tity in question. 

Additionally, care must be taken to avoid feeding 
crosswind gusts directly into the lateral control laws 
via sideslip feedback. High-speed flight through at- 
mospheric disturbances could otherwise cause sharp 
peaks in commanded surface deflections. Since ex- 
tremely sharp peaks in the commanded position of 
control surfaces are generally undesirable and tend 
to cause actuator-rate saturation, a low-pass filter is 
added to the sideslip-angle feedback loop. The filter 
is complemented with derivative feedback to elimi- 
nate its effect on the stability of the closed loop and 
to prevent it from interfering with the response to 
pilot or guidance system commands. Because of the 
derivative feedback through the filter, this network is 
referred to as a “complementary filter.” The comple- 
mentary filter has a unity transfer function and does 
not affect the stability of the closed-loop system. The 
filter acts only to reduce high-frequency control ac- 
tion in response to atmospheric disturbance inputs, 
and not to commanded control inputs. It would also 
be necessary to filter the signal appropriately to elim- 
inate corruption of the measurement as a result of 
structural vibrations if a flexible vehicle model were 
used. A formulation of the complementary filter is 
included in appendix A. Several factors are also in- 
cluded in this appendix that reflect important consid- 


erations in inner-loop controls design for hypersonic 
vehicles. Although a relatively simple inner-loop de- 
sign technique has been used for this investigation, 
more sophisticated designs can also be used in con- 
junction with the resolver concept presented in this 
report. 

Outer-Loop Flight Controls 

The outer loops of the control system shown in fig- 
ure 2 track guidance commands in altitude and cross 
range. They convert vertical and lateral position er- 
rors into vertical and lateral acceleration commands. 
For the outer loops, a proportional-integral structure 
with derivative feedback is used to provide damping 
(p-i-d control). A dual-loop architecture that reg- 
ulates altitude and cross range is used to provide 
coordinated tracking of maneuver commands from 
the guidance system. The outer loops are synthe- 
sized separately, but their command signals to the 
inner-loop controls are coupled through the action 
of the resolver presented subsequently in this report. 
Inner-loop dynamics are assumed to be sufficiently 
fast so as not to invalidate the integrated system de- 
sign. This is not an unreasonable assumption, since 
trajectory parameters are slow to respond in hyper- 
sonic flight. Figure 3 is a block diagram of the basic 
altitude loop architecture. The feedback gains in the 
p-i-d control loop may be specified in terms of fre- 
quency and damping as described below. 

The inner loop for vertical acceleration response is 
assumed to operate with sufficient speed to be trans- 
parent to the outer-loop dynamics, so its transfer 
function is represented as a gain of 1. The closed- 
loop transfer function for the altitude loop shown in 
figure 3 may then be expressed as 

h _ K hd K h ( s + K hl) / 2 n 

^cmd T I^hd 8 ^ T A/j^A/jS + A/j^A/jA/jJ 

The third-order characteristic equation of this trans- 
fer function can be represented in terms of a product 
of first- and second-order expressions as 

(s 2 + 2 (h^h s + ( s + °7 1) = 0 (3) 

or 

S 3 + (°7i + 20 I^h) S 2 + ( 20 I^h^h + s + = 0 

( 4 ) 

where and represent frequency and damping of 
the harmonic response, and is the real root asso- 
ciated with the exponential portion of the response. 
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Feedback gains Ky, and Kyj for the p-i-d con- 

troller depicted in figure 3 are now solved in terms 
of frequency and damping of the second-order char- 
acteristic, and the real first-order root, by equating 
coefficients in equation (4) with coefficients in the 
denominator of equation (2) as follows: 


K hd = °7i + 2 ( h Lv h 


“>l + 2Q l w h a h 

a h + 20 i^h 
U h a h 

w h + 2 Ch a h 


(5a) 

(5b) 

(5c) 


The lateral loop has an identical p-i-d architec- 
ture, and the feedback gains are determined by anal- 
ogous equations; the subscript h is replaced by the 
subscript y in equations (2) to (5). This formulation 
allows specification of the outer-loop pole constella- 
tion in terms of frequency, damping, and real root lo- 
cation for the altitude regulation circuit and for the 
cross-range regulation circuit. These p-i-d control 
loops for altitude and cross range are interconnected 
through the resolver as shown in figure 4. They 
produce a vertical and lateral acceleration command 
combination (/i cmd , Ocmd) that the resolver uses to 
define a normal load factor and bank angle command 
combination (ra cmd , <?Wt)- 

Interface Between Inner and Outer 
Loops 

The resolver acts as an interface between the 
inner- and outer-loop control laws described thus far. 
It couples the altitude and cross-range feedback loops 
by converting vertical and lateral acceleration com- 
mand combinations into an aerodynamic lift vector 
command defined by a normal load factor and bank 
angle. Consider the coordinate system defined in fig- 
ure 5, where n v and n/ are vertical and lateral load 
factors in g units with respect to an Earth reference 
frame. (See ref. If.) For small flight-path angles, 
n v and rif are simply components of the normal load 
factor n, given by n v = ncoscj) and rif = rasing. 
Point A on figure 5 corresponds to straight and level 
flight, since the lateral load factor component n f is 
zero and the vertical load factor component n v is 
equal to that required to maintain constant altitude. 
Any condition along the horizontal line that passes 
through point A denotes a level turn, since n v at 
any point along this line is equal to the vertical load 
factor component required to maintain constant alti- 
tude. The point designated by B, referred to as the 
“ballistic origin” throughout this report, denotes the 


condition at which the vehicle normal load factor n is 
zero (zero lift). The vertical distance that separates 
points A and B is equal to I — V 2 /Rg, where V is 
the vehicle inertial velocity, R is its altitude measured 
from the Earth’s center, and g is gravity. The term 
V 2 / Rg represents an effect referred to as centrifugal 
relief. It results from the centripetal acceleration of 
the vehicle as it flies about a spherical Earth. Lines 
of constant load factor appear as concentric circles 
about the ballistic origin as depicted in figure 5. In 
subsonic flight, the centrifugal relief term is negligi- 
ble and the horizontal line through point A intersects 
the vertical axis at a load factor of lg (i.e. , a vehicle 
in straight and level flight experiences lg). In hy- 
personic flight, the relief term becomes appreciable 
and causes location A to move closer to the ballis- 
tic origin, B. Therefore, the aircraft no longer needs 
to produce a load factor of lg to maintain straight 
and level flight. At extremely high Mach numbers, 
the centrifugal relief term can approach lg. When 
V 2 / Rg = 1, orbital velocity has been achieved, so no 
aerodynamic lift is required to maintain level flight, 
and point A coincides with point B in figure 5. 

It is possible to define a rectangular-to-polar 
coordinate transformation that converts desired ver- 
tical and lateral accelerations (/i cmd , g cmd ) into a 
normal load factor and bank angle command com- 
bination (ra cmd , </> cmd ); this combination represents 
polar coordinates about the ballistic origin. The 
commanded vertical and lateral accelerations are for- 
mulated based on altitude and cross-range errors, as 
previously described in the outer-loop controls sec- 
tion of this report. Desired load factor components 
may be expressed in terms of these vertical and lat- 
eral acceleration commands as follows: 

h 1/2 

w v,cmd — hi - ~j~ (6a) 

a Ra 

2/cmd / c i \ 

n l, cmd = — — (6b) 

Consider a point defined by a given value of 
(m^cmdi ra t,cmd) as depicted in figure 5. The trans- 
formation to polar coordinates is given by 


n = \/ ( n v ) 2 + {rif) 2 (7a) 

and 

4> = arctan ( — ] (7b) 

\n v J 

If the aircraft comes to the angle of attack required 
to achieve the normal load factor n and rolls to the 
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bank angle <f>, then the desired vertical and lateral 
accelerations (h cmc y cm( i) are produced. 

It is apparent from figure 5 and from equa- 
tions (7a) and (7b) that this coordinate transforma- 
tion encounters practical difficulties when (n v , ng) 
command combinations lie extremely near the ballis- 
tic origin. In this vicinity, lift is relatively small, so 
large bank angles are required to effect minor changes 
in cross range. An (n v , ng) command sequence that 
passes directly through the ballistic origin B would 
produce a discontinuity in commanded bank angle of 
180°. For the example design presented subsequently, 
lower limits on allowable angle-of-attack values spec- 
ified by engine tolerances result in a minimum al- 
lowable load factor and thereby avoid the problems 
that would be encountered by operating in extremely 
close proximity to the ballistic origin. In this de- 
velopment, it is assumed that no thrust vectoring is 
available to augment the lift force and that the only 
means of achieving a given («„ cm d, ng cmd) combina- 
tion is by modulating angle of attack to produce the 
desired normal load factor n while rolling the vehicle 
through bank angle <f> to redirect the lift vector. 

Consideration of Flight Condition 
Constraints 

Not all combinations of («„ cm d, ng cmd) produce 
physically achievable or practical values of n and (f>. 
It is therefore necessary to consider boundaries that 
specify a region of allowable («„ cm( i, ng cm d) com- 
binations in the resolver coordinate space. Current 
hypersonic propulsion system concepts are sensitive 
to variations in angle of attack, angle of sideslip, and 
dynamic pressure. Recent studies indicate that dra- 
matic variations in thrust can result from angle-of- 
attack perturbations as large as 2° (ref. 10). For this 
reason, it is assumed that propulsion system toler- 
ances result in an upper and lower limit on angle of 
attack about some nominal operating condition. For 
a given dynamic pressure, the constraints on angle 
of attack result in normal load factor limits, ra max 
and ra m ; n , that would appear as concentric circles 
about the ballistic origin in figure 5. These limits 
are shown in figure 6, which portrays one quadrant 
of the coordinate system from figure 5. The use of 
these load factor limits to define a region of allowable 
( ra u,cmcb n l,c md) combinations in the resolver coordi- 
nate transformation assures that commanded angles 
of attack never violate the constraints, regardless of 
the magnitude of the trajectory error signal. 

If an (n v , ng) command combination lies outside 
the boundary prescribed by ra max in figure 6, it is 
possible to determine the lift and bank angle that 
correspond to this input and then simply limit the 


load factor command to ra max . This limiting scheme 
would result in a radial mapping of all points outside 
the allowable region onto the boundary of a circle of 
radius ra max about the ballistic origin. However, this 
type of limiting scheme produces some undesirable ef- 
fects. For example, a pure cross-range command that 
violates the load factor constraints is represented by 
point 1 in figure 6. (Any point along the dotted line 
in fig. 6 corresponds to a level-turn command, since 
n v is equal to that required to maintain constant al- 
titude.) If the load factor is simply limited to ra max 
while keeping the bank angle constant, the resulting 
coordinates would be represented by point 2 in fig- 
ure 6. Since n v at point 2 is less than that required 
for level flight, the aircraft would descend. Such an 
uncommanded loss of altitude is undesirable, partic- 
ularly if the guidance system is attempting to follow 
some maximum dynamic pressure boundary. There- 
fore, it is clear that simply limiting the load factor 
command issued by the resolver produces an unde- 
sirable effect. 

To maintain the commanded altitude in the pres- 
ence of a limited («„ cm d, ng cmd) value, the points 
outside the allowable region are mapped horizontally 
onto the boundary instead of radially. This mapping 
can be achieved by solving for a maximum allow- 
able lateral component of load factor ng^ max in terms 
of the commanded vertical component of load factor 
ra u,cmd- H the point lies outside the boundary and 
K»,cmd < ™max, then 7i^ max is obtained by solving 
the equation for a circle of radius ra max about the 
ballistic origin B as follows: 


n l , max y “max 2 n v,cmd^ ($ a ) 

A pure cross-range command that violates the 
load factor constraints is represented by point 1 in 
figure 7. By using equation (8a), this command is 
mapped onto the location designated by point 2. It 
is clear that no altitude loss would result from the 
use of this limiting scheme. If the point lies outside 
the boundary, and n v cmc \ > ra max (point 3 in fig. 7), 
then the command is mapped onto the location desig- 
nated by point 4. Therefore, the resolver brings the 
vertical acceleration command to a value less than 
Hmax before pursuing the lateral error. In this way, 
priority has been assigned to the altitude error over 
the cross-range error. Giving priority to altitude reg- 
ulation reduces the amount of activity required to 
compensate for the effect of density variation with 
altitude on dynamic pressure. This limiting strategy 
is referred to herein as “altitude priority.” 

The preceding discussion provides for the inclu- 
sion of an upper load factor limit ra max , which results 
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from the maximum angle-of-attack constraint as a 
result of propulsion system tolerances. However, the 
lower limit on angle of attack also implies the exis- 
tence of a minimum allowable load factor ra m ; n . Some 
(n v , rif ) command combinations lie within the radius 
prescribed by the minimum load factor, such as the 
descent command depicted by point 5 in figure 7. To 
achieve the commanded descent rate without violat- 
ing the minimum load factor constraint, the resolver 
invokes a strategy whereby the vehicle is banked to 
attain the desired vertical component of load factor 
ra u,cmd- Therefore, when the command lies within 
the radius prescribed by ra m ; n , the vertical command 
is pursued at the expense of incurring a transient 
cross-range error. This practice is referred to as a 
“bank-to-dive” strategy throughout this report. The 
minimum lateral component of load factor ra/ )ln ; n re- 
quired to enable the descent is defined by the follow- 
ing expression that is analogous to equation (8a): 



By using equation (8b), the pure descent command 
designated by point 5 in figure 7 is mapped onto the 
bank angle and load factor combination designated 
by point 6. The bank-to-dive strategy is consistent 
with the concept of altitude priority. The use of 
this strategy results in descending “S-turns” in re- 
sponse to large descent commands. A set of condi- 
tional statements used to implement this strategy, 
along with a sketch of their corresponding regions 
in the resolver coordinate system, is presented in 
appendix B. 

The resolver mapping scheme described above 
and depicted in figure 7 transforms vertical and 
lateral acceleration commands from the outer loops 
into a normal load factor and bank angle command 
combination to drive the inner loops. During the 
transformation, load factor limits, altitude-priority, 
and bank-to-dive features are incorporated. 


the nominal path. For this reason, a throttle con- 
troller was designed to cancel out the effect of density 
and drag changes on dynamic pressure while the vehi- 
cle executes maneuvers commanded by the guidance 
system. This controller is derived from an expression 
for dynamic pressure rate that includes the effect of 
density variation with altitude as follows: 

q = p 0 V 0 V + ^p=^(T-D) + ^-C p h (9) 
2 m 2 

The throttle controller uses measurements from a 
longitudinal accelerometer and includes estimates for 
drag and dynamic pressure error. No dynamic pres- 
sure sensor is used. Instead, the dynamic pressure 
error is obtained by subtracting an estimate based 
on velocity and altitude from the nominal dynamic 
pressure. A derivation of the throttle control law is 
provided in appendix C. 

Hypersonic propulsion systems are likely to be 
sensitive to excessive variation about the nominal op- 
erating condition, so large transients in the throttle 
command signal are undesirable. For this reason, 
limits are placed on magnitudes of the thrust per- 
turbations commanded by the throttle control law to 
regulate dynamic pressure. 

Implications of Thrust and 
Angle-of- Attack Limits 

Reference has been made to the existence of lim- 
its on the allowable angle-of-attack and throttle per- 
turbations. The magnitudes of such limits have a 
significant impact on the maneuver envelope of a hy- 
personic vehicle. It is possible to develop expressions 
for maximum climb rate, descent rate, and turn rate 
in terms of these limits. In the following paragraphs, 
these maneuver-rate expressions are presented and 
their implications are discussed. 

First, the limits T max , T min , a max , and a min are 
defined as 


Throttle Control Law 

Hypersonic vehicles may follow trajectories that 
are defined along a flight condition boundary of con- 
stant dynamic pressure. Since high dynamic pres- 
sures are generally desirable from the standpoint of 
propulsive efficiency, and low dynamic pressures are 
more desirable from the standpoint of aerodynamic 
heating and thermal management, the trajectory will 
probably follow a dynamic pressure boundary that 
represents a compromise between these two factors. 
It may also be desirable to maintain constant dy- 
namic pressure during hypersonic maneuvers about 


T max 

E?| 

II 

+ A 

(-) 

1 

W 

W 


\w) 

upper limit 

T ■ 

± mm 

_ To 

- A 

(L) 

1 

w 

W 

\w) 

lower limit 

^max 

= op 

+ Act 

upper 

limit 


a min 1 1 u Acq ower j m j t 


(10a) 

(10b) 

(Ha) 

(lib) 


The limits on angle of attack imply upper and 
lower bounds on load factor, ra max and ra m ; n , which 
can be expressed as functions of the aerodynamic 
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model of the vehicle. Assuming a linear lift curve, 
the load factor limits are 


H max 


ra min 


qS 

W 

qS 
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(12a) 
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The existence of T max , together with the desire 
to maintain constant dynamic pressure, suggests an 
upper limit on allowable climb rate commanded dur- 
ing excursions from cruise flight conditions. There is 
some climb rate beyond which the vehicle cannot ac- 
celerate rapidly enough to compensate for the effect 
of density lapse with altitude to maintain dynamic 
pressure. By determining this maximum climb rate 
and then applying it to the altitude command input 
as a rate limit, climb-rate commands that cause vio- 
lation of the dynamic pressure constraint are avoided. 
Along a constant dynamic pressure ascent (or de- 
scent), climb rate is related to thrust and drag by 
solving the dynamic pressure rate (eq. (9)) for h as 
follows: 


q = — (T -D) + %C p h = 0 (13) 

m 1 

k = ^k( T ~ D ( ( 14 ) 

Since this climb-rate expression is based on the 
ability of the vehicle to accelerate or decelerate, 
drag appears in the expression. To determine the 
maximum and minimum climb rates that can be 
accommodated, it is necessary to develop expressions 
for the maximum and minimum anticipated drag, 
Z?max an d Z? m in > based on the vehicle aerodynamics 
and angle-of-attack limits. Assuming a parabolic 
drag polar, the maximum and minimum anticipated 
drag can be approximated as 

-Dmax = qS (Cfyo + Co,a a max + Co, a a max ! ) (15a) 

Anin = qS (C]J,0 + Co, a ^min + Co, a ^min *) (15b) 

The positive climb-rate command limit is set for the 
condition at which the throttle is at maximum in an 
attempt to compensate for the effect of decreasing 
density on dynamic pressure, while the resolver is 
calling for maximum load factor, which results in 
the maximum angle of attack and maximum drag 
as follows: 
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This limit is somewhat conservative, since greater 
climb rates could be commanded without violating 
the dynamic pressure constraint when drag is not at 
the maximum value. Similarly, the negative climb- 
rate command limit (maximum descent rate) is set 
for the condition at which the throttle is at minimum 
in an attempt to compensate for the effect of increas- 
ing density on dynamic pressure, while the resolver 
is calling for minimum load factor, which results in 
the minimum angle of attack and minimum drag as 
follows: 


~2po 

mVoCp 


(r n 




(16b) 


This limit is also conservative, since greater de- 
scent rates could be commanded without violating 
the dynamic pressure constraint when drag is not 
at the minimum value. Additionally, the load factor 
limit implies an upper limit on the turn rate that can 
be achieved. To avoid commands that are in excess 
of this maximum achievable turn rate, a rate limit is 
applied to the heading command input channel. The 
maximum achievable turn rate at a given flight con- 
dition is designated as ip max . The turn-rate equation 
in terms of load factor, bank angle, and velocity is 

<f> (17) 

The maximum turn-rate command limit is set at the 
maximum load factor and maximum allowable bank 
angle as follows: 

An ax = y Sin ()) max (18) 

This equation expresses a rate command limit that 
is applied to changes in commanded heading angle. 
The maximum allowable bank angle used in equa- 
tion (18) is depicted in figure 8. At larger bank an- 
gles, it would no longer be possible to achieve /i C md 
without exceeding the load factor limit. From the 
figure, it is clear that the maximum allowable bank 
angle for use in equation (18) is given by 


An ax — COS 


A _ A 
-ll 1 Rg 


^cmd 


(19) 


In the control system described in this report, 
the limits on climb-rate and turn-rate commands are 
implemented as shaping blocks in the input signal 
paths for altitude and heading immediately following 
the guidance system block in figure 2. In fact, it is 
probable that the guidance system itself would apply 


10 



such limits rather than simply issuing step commands 
in altitude, cross range, or heading. 

The expressions presented in this section illus- 
trate the manner in which limits on allowable angle- 
of-attack and throttle perturbations, together with 
the desire to maintain constant dynamic pressure, in- 
fluence the maneuver envelope of a hypersonic vehi- 
cle about a steady cruise flight condition. Some ma- 
neuvers, however, may represent transitions between 
steady cruise flight conditions. In such cases, the 
nominal throttle setting and angle of attack might be 
allowed to vary during the maneuver. However, sim- 
ilar limiting expressions may still apply based on the 
levels of perturbations that could be tolerated about 
some smoothly varying nominal. The turn-rate and 
climb-rate limits would then express the acceptable 
levels of departure from the transition flight path. 

Design Example and Test Simulation 

To test the control concept presented in this pa- 
per, an unpiloted simulation was developed by us- 
ing a representative hypersonic single-stage-to-orbit 
configuration as an example vehicle. A sketch of this 
hypersonic winged-cone configuration is presented in 
figure 9. The fuselage has a 5° conical forebody 
with wraparound engine nacelles and a 75° swept 
delta wing. The configuration includes a rudder and 
elevons for vertical and lateral control. Table I is a 
summary of geometric characteristics for this config- 
uration. A considerable aerodynamic data base has 
been generated for this configuration as a result of 
extensive wind-tunnel testing and analytical inves- 
tigations with the aerodynamic preliminary analy- 
sis system (APAS) of reference 12. The geometric, 
propulsive, and aerodynamic data for this configura- 
tion are detailed in references 13 and 14. Aerody- 
namic data are included in tables II and III. The 
inner-loop, outer-loop, and resolver elements were 
designed by using this aerodynamic and propulsion 
model for the example vehicle. 

For the purposes of this study, the configuration 
was trimmed at four different hypersonic flight con- 
ditions along a 2000 lb/ft 2 dynamic pressure trajec- 
tory generated by the program to optimize simulated 
trajectories (POST) discussed in reference 15. Con- 
ditions were selected at Mach numbers of 6, 10, 15, 
and 20. Trim values of angle of attack, velocity, and 
altitude are presented in table IV for these four flight 
conditions. Inertial quantities and trim parameters 
that correspond to these four nominal cruise condi- 
tions were used in the feedback gain and limit expres- 
sions for all control system elements presented pre- 
viously in this report and in appendix A. Frequency 
and damping parameters used in the design example, 


and the equation numbers to which they apply, are 
presented in table V. The control design can be sim- 
ilarly adapted to any desired configuration by sub- 
stituting the appropriate vehicle constants into these 
expressions and selecting the desired frequency and 
damping parameters. 

To implement the climb-rate and turn-rate lim- 
its presented previously, it was necessary to esti- 
mate the magnitude of allowable thrust and angle-of- 
attack perturbations about the trim settings. For the 
example design, angle-of-attack perturbations about 
the trim value commanded to perform the given 
maneuvers were limited to ±0.4°. The magnitudes 
of such limits for an actual working vehicle are a 
topic of current research, and these values may or 
may not be overly restrictive. However, one ob- 
jective of this study is to investigate the manner 
in which propulsion and flight condition constraints 
impact the maneuvering capabilities of hypersonic 
aerospacecraft. These effects may be investigated 
without exact knowledge of the constraint magni- 
tudes, so this simulation has been conducted in the 
interest of revealing some issues of importance to this 
general class of vehicles. For this design example, 
throttle limits were placed at —0.10 and 0.30 about 
the trim thrust-to-weight ratio for straight and level 
cruise. According to the propulsive data base pre- 
sented in reference 13, these variations are well within 
the range achievable by varying the fuel equivalence 
ratio for a representative vehicle. Therefore, the limit 
quantities on thrust and angle-of-attack perturba- 
tions for this example vehicle are defined as 
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An accurate model for the dynamics of a propul- 
sion system to be used on a hypersonic air-breathing 
vehicle was not available, but very rapid response 
was assumed possible for this type of engine. For 
this example, a perfect engine model was assumed 
in which the commanded thrust was equal to the 
achieved thrust. Also, the model did not include 
thrust variations with angle of attack. 

A test simulation of the resolver control system 
and the example vehicle was constructed by using a 
control design and analysis software tool that allows 
simulation of all the nonlinear elements described in 
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this report. In this batch simulation, the resolver was 
used to drive the inner-loop controls directly, instead 
of driving a cockpit flight director as would be the 
case in a piloted simulation; also, perfect angle-of- 
attack control was assumed, in which the commanded 
angle of attack was equal to the achieved angle of 
attack (short-period dynamics were omitted). Time 
histories were produced to test the response of the 
controlled vehicle to various guidance commands and 
atmospheric disturbances. In this way, the ability of 
the controller to complete a specified guidance task 
was evaluated. 

A series of example guidance commands were 
generated to simulate the execution of several ma- 
neuvering tasks. Guidance inputs to the controller 
consisted of vertical offset commands, lateral offset 
commands, and heading changes. Responses to sev- 
eral command combinations, as well as a maneuver 
that utilized the bank-to-dive feature of the resolver, 
were evaluated. Responses to three forms of atmo- 
spheric disturbances were also investigated. The hy- 
personic simulation was subjected to headwind gusts, 
crosswind gusts, and density variations. It is im- 
portant to examine the vehicle response to such at- 
mospheric disturbances, since the disturbances may 
give rise to angle-of-attack and angle-of-sideslip ex- 
cursions, which are a major concern in designing con- 
trollers for hypersonic vehicles. All time histories 
presented in this report represent variations about 
the trim values shown in table IV. Table VI is a sum- 
mary of the response plots in terms of figure number, 
input type, and flight condition. Response time his- 
tories from the batch simulation are discussed in the 
following sections. These sections are intended to 
examine the trends and effects apparent in the re- 
sponse time histories at Mach numbers of 6, 10, 15, 
and 20, and to relate these effects back to their phys- 
ical causes. 

Responses to Guidance Commands 

Simulated perturbation responses to an altitude 
step command of 2000 ft for Mach numbers of 6, 
10, 15, and 20 are shown in figure 10. The alti- 
tude responses reveal a large variation in the times 
required to achieve the commanded altitude at the 
different flight conditions. The length of time re- 
quired to achieve the altitude command is primarily 
a function of the climb-rate limits presented in equa- 
tions (16). The variation in speed of response at the 
four Mach numbers is mainly caused by the differ- 
ing values of maximum thrust (T max ) at each flight 
condition and is therefore a result of the manner in 
which limits were implemented in this specific exam- 
ple design. As a result of observing the climb-rate 


command limit, dynamic pressure variations during 
the maneuver were maintained to ±2 lb/ft 2 about the 
nominal 2000 lb/ft 2 condition. The angle-of-attack 
responses presented in figure 10 indicate an initial in- 
crease above the trim value to initiate the climb, fol- 
lowed by a marked decrease to a new trim value at the 
higher altitude. This trim angle-of-attack decrease in 
response to the altitude change along a constant dy- 
namic pressure trajectory was caused by two factors 
of particular significance in hypersonic flight. First, 
the density lapse with altitude required the vehicle 
to accelerate to maintain constant dynamic pressure; 
this effect was apparent in the velocity time histories 
of figure 10. Second, it was stated previously that the 
effect of centrifugal relief in hypersonic flight is to re- 
duce the load factor required to maintain constant 
altitude as velocity increases (eq. (1)) as follows: 


Since the vehicle increased velocity to maintain 
dynamic pressure at the higher altitude, centrifugal 
relief was greater at the new trim condition. There- 
fore, the new trim condition required less lift and a 
lower angle of attack to maintain straight and level 
flight. 

The throttle variations required to regulate dy- 
namic pressure during the maneuver are presented 
in terms of thrust-to-weight ratio in figure 10. The 
Mach 15 and Mach 20 cases reached the maximum 
thrust-perturbation limit during the ascent. Trim 
throttle settings changed slightly in response to the 
lower angle of attack at the end of the maneuver. 

Perturbation responses to a cross-range step com- 
mand of 20 000 ft are shown in figure 11. The cross- 
range time histories are identical for all Mach num- 
bers, since the outer loop was designed to the same 
frequency and damping in each case. Although cross- 
range time histories are identical for the four flight 
conditions, their corresponding ground tracks differ 
dramatically since the trim velocities are dissimilar. 
(See table IV.) The responses in figure 11 show that 
much larger bank angles were required to produce a 
given lateral acceleration at higher velocities; these 
responses illustrate the influence of centrifugal relief 
on high-speed turns. The heading responses in fig- 
ure 11 show that, as expected, a given lateral accel- 
eration produced a larger change in heading at lower 
velocities. As the /3 response indicates, the inner-loop 
controls were designed to keep sideslip to a minimum 
during lateral maneuvers. Altitude variations were 
kept to within 10 ft of the trim setting, and angle- 
of-attack variations required to produce the desired 
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load factors commanded by the resolver did not reach 
the 0.4° limit. Small thrust and velocity variations 
were produced by the throttle control law as it acted 
to regulate dynamic pressure during the maneuver. 

A polar plot of load factor versus bank angle in 
the resolver coordinate system for the 20 000-ft cross- 
range offset maneuver is shown in figure 12. It is 
clear from figure 5 that the time histories in figure 12 
represent excursions of bank angle and load factor 
commands along the horizontal line for level turns 
at each of the four flight conditions. The dramatic 
reduction in vertical load factor required to maintain 
level flight at the higher Mach numbers is apparent. 
This reduction in vertical load factor as a result of 
centrifugal relief causes an increase in bank angle 
required to produce a given lateral acceleration in 
level flight. (See fig. 11.) 

Responses to a vertical command of 2000 ft com- 
bined with a lateral offset command of 20 000 ft are 
shown in figure 13. The resolver produced what is 
essentially a superposition of the responses from fig- 
ures 10 and 11. Dynamic pressure was regulated 
to within ±2 lb/ft 2 of nominal during this maneu- 
ver. All the trends that were apparent in figures 10 
and 11 as a result of the flight condition constraints 
and centrifugal relief effects are also apparent in these 
responses. 

Commanded changes in heading angle were rep- 
resented by ramp inputs in the cross-range command 
signal path. The maximum turn-rate command limit 
presented in equation (18) was applied to this input 
signal. Responses to a commanded heading change of 
10° are shown in figure 14. The variation in speed of 
response is a result of the effect of different trim ve- 
locities on the maximum turn-rate command in equa- 
tion (18). Exceeding this turn-rate command limit 
would have resulted in an uncommanded loss of alti- 
tude in response to the 10° heading change. Again, 
the bank angle required to produce a given lateral 
acceleration increased with Mach number. Altitude 
perturbations were kept to a minimum, and angle-of- 
attack variations reached the upper limit of 0.4° for 
those flight conditions that required sustained oper- 
ation at the maximum allowable turn rate. Clearly, 
at the higher Mach numbers, even minor heading 
changes can require large amounts of time. The 
throttle controller produced minor thrust variations, 
which compensated for the higher drag at larger an- 
gles of attack, to regulate dynamic pressure through- 
out this heading-change maneuver. 

The concept of a bank-to-dive maneuver was de- 
scribed previously in this report. A bank-to-dive ma- 
neuver allows the vehicle to remain at or above the 


minimum allowable angle of attack while achieving 
a high rate of descent. By rolling the lift vector off 
of vertical and maintaining the minimum allowable 
angle of attack, a faster rate of descent is realized 
than would otherwise be possible. Because of the 
significant effect of density variation with altitude, 
it was necessary to relax the constraint of maintain- 
ing constant dynamic pressure to execute a bank-to- 
dive maneuver for the hypersonic vehicle used in this 
example. Since density increases as the vehicle de- 
scends, it is necessary to reduce speed to maintain 
constant dynamic pressure. The maximum allowable 
descent-rate command presented in equation (16b) 
was derived to prevent descent commands that would 
require the vehicle to decelerate faster than possible. 
However, to enter the region of figure 7 where the 
commanded load factor was less than ra m ; n , it was 
necessary to allow a descent-rate command in excess 
of the value expressed in equation (16b). Time histo- 
ries for a descent command step of 5000 ft at Mach 15 
with the bank-to-dive feature are shown in figure 15. 
This response was compared with a response to the 
same descent command without relaxation of the dy- 
namic pressure constraint and consequential descent- 
rate limit. When the descent-rate command limit 
was not relaxed, the bank-to-dive feature of the re- 
solver was never invoked, so a pure vertical descent 
resulted (no lateral activity). It is apparent from 
figure 15 that the bank-to-dive maneuver produces 
a much faster response at the expense of incurring a 
transient cross- range excursion and dynamic pressure 
error. The vehicle was unable to decelerate rapidly 
enough to remain on the dynamic pressure bound- 
ary during the bank-to-dive descent. The cross-range 
and bank angle excursions were plotted on a time 
scale of 200 sec rather than 800 sec, since no lat- 
eral activity occurred for the remainder of the ma- 
neuver. The lateral flight-path deviation remained 
within ±3000 ft, and the bank angle excursion was 
about ±50°. The trim angle of attack at the end of 
the maneuver was higher than at the start of the ma- 
neuver. The deceleration reduced the centrifugal re- 
lief term, V 2 /Rg, so the vehicle had to produce more 
lift to remain in level flight. (The dynamic pressure 
was the same at the beginning and end of the ma- 
neuver, so it is not a matter of requiring a greater 
angle of attack to produce the same amount of lift.) 

A comparison of the bank-to-dive responses to the 
5000-ft descent command for Mach 6, 10, 15, and 20 
is shown in figure 16. Although the altitude per- 
turbations are almost identical, the dynamic pres- 
sure excursions are greater at the faster and higher 
flight conditions. It is possible that the addition 
of a deployable drag-producing device would enable 
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the vehicle to decelerate rapidly enough to utilize a 
bank-to-dive maneuver scheme while satisfying the 
constant dynamic pressure constraint, although the 
need for such a device would be configuration and 
flight condition dependent. Larger bank angle excur- 
sions were again experienced at the faster flight con- 
ditions as a result of greater centrifugal relief. The 
bank-to-dive method might be used during reentry, 
abort, or nonorbital missions. 

Responses to Atmospheric Disturbances 

Responses to three types of atmospheric distur- 
bances were simulated at each of the four hypersonic 
flight conditions. The first type of disturbance was 
a 30-knot headwind with a fO-mile ramp onset. The 
second type of disturbance was a 30-knot crosswind 
with a fO-mile ramp onset. After reaching a peak 
value of 30 knots, the headwind or crosswind re- 
mained constant through the rest of the simulation. 
The third type of disturbance was an atmospheric- 
density perturbation represented as a fO-percent de- 
crease below nominal in the form of one cycle of a 
negative (f — cosine) wave. This disturbance was in- 
tended to represent flight through a region of rarefied 
air in the upper atmosphere. Such density variations 
have been recorded during shuttle reentry flights at 
very high altitudes (ref. 16), although much uncer- 
tainty still exists as to the structure and magnitude 
of density perturbations that a hypersonic aerospace- 
craft might encounter. The negative 10-percent den- 
sity pulse was defined with a horizontal length of 
10 miles. Since these disturbances were defined spa- 
tially (in terms of miles), the vehicle experienced the 
perturbations at different rates in the time domain, 
depending on trim velocity of the flight condition. 
Recent studies indicate that the severity of these dis- 
turbances is light compared with the strength of per- 
turbations that may be encountered in the actual 
environment (ref. 16). Therefore, these responses 
are used to illustrate the general manner in which a 
hypersonic vehicle equipped with the control design 
presented in this report would react to such distur- 
bances; these responses should not be interpreted as 
an exhaustive study of robustness to environmental 
perturbations. 

Responses to the 30-knot headwind with a 10- 
mile ramp onset at each of the four flight conditions 
are shown in figure 17. Very small angle-of- attack 
and altitude variations were experienced. Dynamic 
pressure and thrust variations were plotted on a 
shorter time scale, since these transients occurred 
more rapidly. The throttle control law reduced the 
thrust level to the lower limit to decelerate the vehicle 
as rapidly as possible; this deceleration compensated 


for the effect of the headwind on dynamic pressure. 
As previously mentioned, a model of the engine 
dynamics for the hypersonic propulsion system was 
not readily available for use in this simulation. The 
throttle transients in response to the headwind were 
quite rapid, and an accurate model of the engine 
dynamics might influence this response and possibly 
require some modification or tuning of the throttle 
control law. 

Cross-range deviation and sideslip responses to 
the 30-knot crosswind with a fO-mile ramp onset are 
shown in figure 18. Perturbations in sideslip angle 
were less than a tenth of a degree, and there were 
negligible changes in angle of attack, thrust, bank 
angle, and heading. The resulting cross-range de- 
viations were negligible, although they took longer 
to eliminate. (Note the longer time scale on cross 
range.) At hypersonic speeds, the longitudinal com- 
ponent of velocity greatly diminished the effect of any 
crosswind component. 

Responses to a negative 10-percent density pulse 
with a horizontal length of 10 miles are shown 
in figure 19. Transients that resulted from the 
density pulse were greatest at the slowest flight 
condition. Angle-of-attack and throttle variations 
remained within the perturbation limits specified 
previously in this report for the design example. The 
effect on dynamic pressure was plotted on a time 
scale of 20 sec, since the vehicle traversed the dis- 
turbance so rapidly. It is clear that the throttle 
controller could not act quickly enough to neutralize 
the dynamic pressure excursion produced by such a 
disturbance. 

The spatial dimension of the density pulse was 
varied at the Mach 10 flight condition to examine 
the effect of disturbance duration on the vehicle re- 
sponse. Responses at the Mach 10 flight condition to 
a negative 10-percent density pulse with lengths of 5, 
10, 15, and 20 miles are compared in figure 20. The 
longer pulses produced the greatest transients, al- 
though perturbations again remained within the lim- 
its specified previously in this report. The throttle 
control was again unable to respond with sufficient 
speed to neutralize the effect of density variations 
on dynamic pressure. In the resolver control sys- 
tem, altitude is regulated as a trajectory variable, 
while dynamic pressure is handled by an inner-loop 
throttle control law. This type of control scheme is 
well-suited to regulation of dynamic pressure during 
a planned maneuver, but it cannot provide perfect 
regulation of dynamic pressure in response to atmo- 
spheric density perturbations. To completely cancel 
the effect of density variations on dynamic pressure, 
the vehicle would be required to change velocity or 
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altitude with a faster response than appears to be 
reasonable. The inability to completely compensate 
for the effects of density variations on dynamic pres- 
sure suggests that a practical hypersonic propulsion 
system must be able to tolerate a certain level of dy- 
namic pressure fluctuation. This finding is consistent 
with previous studies that examined various impli- 
cations of regulating different trajectory parameter 
combinations. (See refs. 2 to 6.) 

Concluding Remarks 

A control design option is presented for execut- 
ing coordinated maneuvers in hypersonic flight while 
regulating key parameters such as angle of attack, 
angle of sideslip, and dynamic pressure. The design 
involves the use of a coordinate transformation, re- 
ferred to as a “resolver,” which converts vertical and 
lateral acceleration commands into an aerodynamic 
lift vector command specified by a normal load factor 
and bank angle. The vertical and lateral acceleration 
commands are based on feedback errors from com- 
manded values of altitude and cross range. The con- 
trol system was applied to an example configuration 
and implemented in an unpiloted digital simulation 
at Mach numbers of 6, fO, 15, and 20. The simula- 
tion was used to illustrate the manner in which the 
control system responds to various commands and 
atmospheric disturbances. 

A vehicle in hypersonic flight may approach or- 
bital velocities while still in the atmosphere. The 
centripetal acceleration of the vehicle as it circles the 
Earth then becomes significant enough to cause a 
noticeable reduction in the load factor required to 
maintain constant altitude. This effect, referred to as 
centrifugal relief, was demonstrated in the responses 
of the simulation to commanded altitude, heading, 
and cross-range changes. A dramatic reduction in 
vertical load factor required to maintain level flight 
at the higher Mach numbers was observed. This re- 
duction in vertical load factor caused the bank angle 
that was required to produce a given lateral accel- 
eration to increase with Mach number. Centrifugal 
relief also caused changes in trim angle of attack in 
response to altitude variations along a constant dy- 
namic pressure trajectory. The control design pre- 
sented in this report accounts for the effect of cen- 
trifugal relief while tracking altitude and cross-range 
commands. It also incorporates load factor limits, 
prioritizes altitude regulation over cross-range regu- 
lation, and includes a feature whereby the lift vec- 


tor can be rolled off of vertical to achieve a desired 
descent rate at the expense of incurring a transient 
cross-range error (“bank-to-dive”). 

Several implications of propulsion and flight 
condition constraints for the maneuvering of air- 
breathing hypersonic vehicles were examined. Angle- 
of-attack and throttle perturbation constraints, com- 
bined with centrifugal relief effects and the 
desire to maintain constant dynamic pressure, sig- 
nificantly impact the maneuver envelope for such 
vehicles. Turn-rate, climb-rate, and descent-rate lim- 
its with respect to maneuvering about some hyper- 
sonic cruise condition were expressed in terms of 
these constraints. Density variation with altitude 
strongly influences climb-rate and descent-rate limits 
and throttle modulation if dynamic pressure is to be 
maintained during vertical transitions between cruise 
flight conditions. In particular, the high descent rates 
that are achievable by using a bank-to-dive maneu- 
ver scheme required a departure from the dynamic 
pressure constraint, since the vehicle was unable to 
decelerate rapidly enough to compensate for the ef- 
fect of density lapse on dynamic pressure. How- 
ever, the bank-to-dive method might be used during 
reentry, abort, or nonorbital missions. 

Of the atmospheric disturbances investigated, 
density pulses caused the greatest transient response. 
The findings of this study indicate that hypersonic 
propulsion systems should be designed with the ca- 
pacity to tolerate fluctuations in dynamic pressure 
that result from high-speed flight through atmo- 
spheric density perturbations. The effects of atmo- 
spheric disturbances on advanced propulsion systems 
must be investigated and provided for in any practi- 
cal hypersonic design. 

This research has defined an automatic control 
design for executing coordinated maneuvers in hy- 
personic flight while regulating key flight condition 
parameters. The primary element of this control de- 
sign provides an environment for the incorporation of 
various constraints associated with high-speed flight. 
Additionally, this work has illustrated the manner 
in which propulsion and flight condition constraints 
impact the maneuvering capabilities of air-breathing 
hypersonic aerospacecraft. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
December 11, 1991 
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Appendix A 

Development of Inner-Loop Flight 
Controls 

This appendix describes the derivation of the lat- 
eral inner loops, which were designed to stabilize 
the vehicle dynamics, track bank angle commands 
from the resolver, and maintain zero sideslip dur- 
ing the maneuver. A simple pole-placement tech- 
nique is used that allows specification of the inner- 
loop responses in terms of desired frequency and 
damping. Two separate control loops are synthe- 
sized to drive the rudder and the ailerons. The esti- 
mation of sideslip angle with a lateral accelerometer 
and the concept of a complementary filter to reduce 
high-frequency control action that results from at- 
mospheric disturbances are also described. Although 
a relatively simple inner-loop design technique has 
been used for this investigation, more sophisticated 
designs can also be used in conjunction with the re- 
solver concept presented in this report. 

Rudder Control Law 


force due aileron deflection yields 

sr(s) - Np/3(s) - Ng r [Krp/3(s) + Kr r r(s)\ = 0 

(A6) 


s(3(s) + r(s)-Yp(3(s)-Yg r [Krp(3(s) + Kr r r(s)] =0 

(A7) 


This system can be expressed in matrix form as 


' s — N$ r Kr r —Np — Nf r Krp 


-r(s)- 


'O' 

. 1 - Y Sr Kr r s - Yp - Y Sr Krp . 


J(s)_ 


. 0 . 


This matrix yields the following second-order char- 
acteristic equation for the closed-loop system: 

s 2 - (Yp + Y 6r Krp + N 6 r Kr r ) s 

+ [N f 3 + N ir Kr p + Kr r ( N Sr Yp - Y Sr Np)] = 0 (A9) 

The coefficients of this characteristic equation are ex- 
pressed below in terms of frequency co r and damping 
( r for the second-order system: 


The rudder is used to maintain zero sideslip dur- 
ing a maneuver. First, a Dutch Roll approximation 
is applied to address the yaw control, and only the 
yaw-rate and sideslip equations are given as follows: 

r = N af3 + Npp + N r r + N^ a 6 a + N^ r 6 r — - — p (AI) 

\i = Ypf} + Y p p + Y r r + Y Sa S a + Y Sr S r - r + y cos 0 O sin <j> (A2) 

The final term in equation (A2) involves gravity 
and is negligible in hypersonic flight since velocity is 
very large. Also, the effect of rotary derivatives is of- 
ten small at hypersonic speeds, since velocity appears 
in the denominator of the dimensionalizing expres- 
sion. If the rotary derivative and inertial coupling 
terms in equations (AI) and (A2) are neglected, a 
Laplace transformation yields the following dynamic 
system: 

sr(s) = Np(3(s) + N Sa 6 a (s) + N Sr 6 r (s) (A3) 

sP{s) = Y r J{s) + Y s J a {s) + YgMs) - r{s) (A4) 

A rudder deflection command is defined using yaw 
rate and angle-of-sideslip feedback of the following 
form: 

6 r (s) = Krpf3(s) + Kr r r(s) (A5) 

Substituting the rudder deflection command into ex- 
pressions (A3) and (A4) and neglecting yaw and side 


2 Cr^r = - (Yp + Y Sr Krp + N Sr Kr r ) (AI 0) 

w r = N [3 + N 8r Kr f3 + Kr r (Ng r Yp - Yg r Np) (All) 

Neglecting the relatively small terms, Kr r ( Ng r Yp 
— Yg r Np), and solving for the feedback gains, Krp 
and Kr r , in terms of frequency and damping, yields 


Krp = 


u 2 r ~ Np 
Ng r 


(A12a) 


and 


Kr r 


2 Cr^r + Yp + Yg r R rp 
-Ng r 


(A12b) 


It is now possible to quickly determine the feed- 
back gains required to produce a given closed-loop 
frequency and damping combination. This combi- 
nation allows the designer to place the closed-loop 
poles at a location that might be desirable from a 
handling-qualities perspective. 


Angle-of-sideslip estimation. The reliabil- 
ity of air-data measurements during some phases of 
hypersonic flight may be uncertain. The control 
system designed in this report utilizes an angle-of- 
sideslip feedback that is constructed from a lateral- 
accelerometer reading rather than an air-data sensor 
array. The lateral accelerometer is placed at the in- 
stantaneous center of rotation in response to a rudder 
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deflection; this placement eliminates the need to cor- 
rect the accelerometer reading for rudder effects. It 
would actually also be necessary to filter the signal 
appropriately to eliminate corruption of the measure- 
ment as a result of structural vibrations. The equa- 
tion for the lateral acceleration measured by an ac- 
celerometer located a distance A forward of the center 
of gravity is given by 

Ay, acc = Vo (YpP + Yg r S r ) + A (Npf3 + Ng r 6 r ) 

(A13) 

Solving this expression for /3 yields an angle-of- 
sideslip measurement, in terms of the accelerometer 
reading and the rudder deflection, of 

acc Sr (VoYg r + \Ng r ) Ml-ll 

P ~ V 0 Yp + A Np 1 j 


is added to the sideslip-angle feedback loop. The fil- 
ter is complemented with derivative feedback to pro- 
duce a unity transfer function, which eliminates its 
effect on the stability of the closed-loop system and 
prevents it from interfering with the response to pi- 
lot or guidance system commands. Because of the 
derivative feedback through the filter, this network 
is referred to as a complementary filter. The expres- 
sion for the filtered angle-of-sideslip feedback is 


P 


r/3 + P 

TS + 1 TS + 1 


(A16) 


where /3 is obtained from equation (A14) and /3 is 
given by 

P = Vp(3 + Yg a 6 a + Yg r 8 r — r (A17) 


By placing the accelerometer at the lateral center 
of rotation, it is possible to eliminate the influence 
of rudder deflection on the angle-of-sideslip measure- 
ment. If 


then 


_ -VoY Sr 

N Sr 


(A15a) 


P 


Ay, acc 

VoYp + A N p 


(A15b) 


This expression for sideslip can be substituted into 
equation (A5) to make the lateral inner-loop con- 
troller independent of air-data measurements. The 
expression for A depends on parameters that change 
with flight condition; the lateral center of rotation 
is not constant over the range of hypersonic flight 
conditions experienced by the vehicle. Since it is 
unreasonable to relocate the accelerometer as flight 
condition varies, it is probable that the instrument 
would be placed at some location that represents the 
average position of the center of rotation for those 
flight conditions at which the air-data system mea- 
surements are the least reliable. For the simulated 
time histories presented in this report, which include 
flight conditions at Mach 6, 10, 15, and 20, the lateral 
accelerometer was placed at the center of rotation for 
the Mach 6 flight condition. 


Rotary derivatives have been neglected in equa- 
tion (A17). The time constant can be selected to 
provide a reasonable amount of disturbance rejec- 
tion and a favorable response characteristic. Expres- 
sion (A16) can be substituted for sideslip feedback in 
the rudder control law (eq. (A5)). Ideally, the com- 
plementary filter has a unity transfer function and 
does not affect the stability of the closed-loop system. 
In practice, the degree to which the transfer function 
of equation (A16) differs from unity depends on the 
accuracy of the estimated derivative signal used in 
this equation. To this extent, the complementary fil- 
ter acts to reduce high-frequency control action in 
response to atmospheric disturbance inputs but not 
to commanded control inputs. 

Aileron Control Law 

The roll-control law drives the ailerons to achieve 
bank angle commands issued by the resolver. The 
roll-rate and bank angle equations are as follows: 

P = Lp/3 + L p p + L r r + Lg a 6 a + Lg r 6 r + r (A18) 

J-X 

4> = p + r tan #o cos <t> (A19) 

Neglecting the rotary derivative and inertial coupling 
terms in equation (A18) yields 


Complementary filter. Sudden crosswind 
gusts are fed directly into the lateral control-surface 
deflections via sideslip feedback under the current 
rudder control law. To avoid sharp peaks in the com- 
manded position of control surfaces, a low-pass filter 


p — Lpj3 + Lg a 6 a + Lg r 6 r (A20) 

The simplifying assumption is made that <f> « p 
in equation (A19). A feedback control law for 
aileron deflection that allows decoupling of yawing 
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and rolling motions by eliminating the effect of rud- 
der deflection and angle of sideslip on equation (A20) 
is defined as follows: 

6 a = Kap/3 + Ka Sr 6 r + Ka^ eTT + Ka^<f> (A21) 

In equation (A21), 4> eTT is defined as the commanded 
roll angle minus the actual roll angle. This control 
law also includes roll-attitude and roll-rate feedback, 
which provides a means of specifying frequency and 
damping in roll response. The angle-of-sideslip sig- 
nal used in this expression is the complementary fil- 
tered estimate in equation (A16) from the lateral ac- 
celerometer. The rudder deflection command is also 
used in this control law. Substituting the aileron con- 
trol law into equation (A20) yields 

^ = Lpj3 + L$ r 6 r + Lj a + Ka^ r 6 r + 

(A22) 

To achieve the desired decoupling of roll and yaw 
responses, let 

—Lp. 

I\ a = 3 ! A 23a ! 

LSa 

and 

Ka Sr = (A23b) 

L 8a 


which, upon substitution into equation (A22), yields 

4> = LgaKa^err + L Sa Ka^it> (A24) 

A Laplace transform of equation (A24) yields the 
following second-order characteristic equation: 

s 2 - L Sa Ka ^s - L Sa K(i 0 = 0 (A25) 

The feedback gains in this expression are easily solved 
for in terms of desired frequency and damping of the 
closed-loop system to yield 

2 

Kcla = — j— (A26a) 

~ L Sa 

and 

= ^>a (A2gb) 

^Sa 

These gains, along with those in equations (A23), 
are used in equation (A21) to produce the final 
aileron control law. The rotary terms neglected in 
the derivation of this control law were of little signif- 
icance at the hypersonic flight conditions that were 
simulated. The inherent aerodynamic roll damping 
of the vehicle is negligible at these flight conditions, 
so most of the damping is supplied by the roll-rate 
feedback loop in the aileron control law. 
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Appendix B 

Conditional Statements Used in 
Bank-to-Dive Strategy 

This appendix describes some conditional state- 
ments that are needed to successfully implement the 
bank-to-dive limits in the resolver coordinate system. 
Three regions in the resolver coordinate system are 
illustrated in sketch A. Region 1 consists of all points 
for which the normal load factor is less than ra m ; n , the 
minimum allowable load factor prescribed by equa- 
tion (12b). Region 2 consists of all points outside of 
region 1 for which n„ cmc i is less than ra m ; n . Region 3 
consists of all points outside of regions 1 and 2. Once 


a bank-to-dive maneuver has been initiated by a com- 
mand lying within region 1, a lateral error is incurred 
to achieve the desired descent rate. It is then neces- 
sary to suspend the normal operation of the resolver 
until the acceleration command no longer lies within 
the areas designated as regions 1 and 2 in sketch A. 
Once the acceleration command sequence has passed 
into region 3, the lateral error can be eliminated. 
Omission of these conditional statements would re- 
sult in a bank angle command chatter between left 
and right halves of the resolver coordinate plane each 
time the acceleration command sequence passes from 
region 1 to region 2. In the sketch, B-T-D is a con- 
dition that can be either true or false, depending on 
the location of the current acceleration command. 


Sketch A 
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Appendix C 

Development of Throttle Control Law 

This appendix describes the derivation of a feed- 
back control law for the throttle that has been de- 
signed to regulate dynamic pressure during a ma- 
neuver. The throttle control law is derived from the 
expression for dynamic pressure rate, which is given 

by 


V 2 

q = PoV 0 V + -^-p (Cf a) 

where velocity variations and density lapse with alti- 
tude may be approximated as 


and 


V = 


T-D 

m 


(Clb) 


i,= %-h = C f h 


(Clc 


The term C p is the density lapse coefficient for a given 
flight condition. Substitution of equations (Cfb) and 
(Cfc) into equation (Cfa) yields 

q = ^-T - ^-D + ^-C p h (C2) 
m m 2 


A throttle control law is defined, based on equa- 
tion (C2), to compensate for the effect of altitude 
and drag variations while providing first-order regu- 
lation of dynamic pressure as follows: 

T cmd = D- + Kfq eTT (C3) 

Z P o 


Substitution of T cm( j for T in equation (C2) produces 


v _ PoVq - 
q — -A-qQeTT 
m 


(C4) 


and 

q( s )( s -£M K \ Q (C5) 

\ m J 

which yields the following first-order characteristic 
equation: 


where 


s T (Tq = 0 


-do fo 


m 


(C6a) 

(C6b) 


The dynamic pressure feedback gain can be expressed 
in terms of the real root in the following first-order 
characteristic: 


K 




— mCTq 

do bo 


(C7) 


The throttle controller described in equation (C3) 
includes a drag estimate and a dynamic pressure er- 
ror feedback. Drag is estimated by solving equa- 
tion (Cfb) as follows: 


D = T prev - mVacc (C8) 

where the inertial velocity variation is obtained from 
a longitudinal accelerometer and T prev represents the 
thrust command from the previous time step. The 
dynamic pressure error is obtained by subtracting 
an estimate based on velocity and altitude from the 
nominal dynamic pressure. Therefore, 


?err = ?0 “ <Z (C9) 

where the dynamic pressure estimate is given by 
^ W 2 V 2 r 

q = -y- = — [do + Cp(h - do)] (CfO) 

Substituting the drag and dynamic pressure esti- 
mates from equations (C8), (C9), and (CfO) into 
equation (C3) yields the final throttle control law as 
follows: 


mV(jC„ ■ 

Tcmd = Tprev “ ™C acc - P h + Kg 


qo - -y [do + Cp(h - do)] [> (Cff) 


The real root in equations (C6) should be cho- 
sen to prevent excessive dynamic pressure excur- 
sions without producing extreme throttle transients. 
This development of the throttle control law depends 
heavily on an accurate model of atmospheric density 
variation with altitude, and its performance may be 
strongly influenced by model uncertainty. 

The control format described in this report as- 
sumes that the nominal throttle setting is dictated 
by guidance and trajectory considerations and that 
some perturbations about this nominal setting are 
permitted so that dynamic pressure can be regulated. 
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It is possible to envision formats in which pertur- 
bations about the nominal throttle setting are not 
permitted. In such a case, the guidance system may 
command altitude changes to regulate dynamic pres- 


sure in accordance with the throttle setting. The 
resolver concept presented in this report may still 
be used in conjunction with such alternative control 
formats. 
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Table I. Geometric Characteristics of Winged-Cone Configuration 

Wing: 

Reference area (includes area projected to fuselage centerline), ft 2 3603.0 

Aspect ratio f.O 

Span, ft 60.0 

Leading-edge sweep angle, deg 76.0 

Trailing-edge sweep angle, deg 0.0 

Mean aerodynamic chord, ft 80.0 

Airfoil section Diamond 

Airfoil thickness-to-chord ratio, percent 4.0 

Incidence angle, deg 0.0 

Dihedral angle, deg 0.0 

Wing flap (elevon): 

Area, each, ft 2 92.3 

Chord (constant), ft 7.2 

Inboard-section span location, ft 15.0 

Outboard-section span location, ft 27.8 

Vertical tail: 

Exposed area, ft 2 645.7 

Theoretical area, ft 2 1248.8 

Span, ft 32.5 

Leading-edge sweep angle, deg 70.0 

Trailing-edge sweep angle, deg 38.1 

Airfoil section Diamond 

Airfoil thickness-to-chord ratio, percent 4.0 

Rudder: 

Area, ft 2 161.4 

Span, ft 22.8 

Ratio of chord to vertical tail, percent 25.0 

Axisymmetric fuselage: 

Theoretical length, ft 200.0 

Cone half-angle, deg 5.0 

Cylinder radius (maximum), ft 12.9 

Cylinder length, ft 12.9 

Boattail half-angle, deg 9.0 

Boattail length, ft 40.0 

Moment reference center, ft 124.0 
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Table II. Aerodynamic Stability and Control Derivatives for Design Example 



Mach 6 

Mach 

10 

Mach 

15 

Mach 20 

Cyp 

-9.34 

X 

10“ 3 

-7.34 

X 

10“ 3 

-6.52 

X 

10“ 3 

-6.20 

X 

10“ 3 

Cy Sr 

2.40 

X 

10“ 4 

1.30 

X 

10“ 4 

9.00 

X 

10“ 5 

7.00 

X 

10“ 5 

c b 

-7.00 

X 

10“ 4 

-4.70 

X 

10“ 4 

-3.40 

X 

10“ 4 

-2.70 

X 

10“ 4 

C kr 

7.00 

X 

10“ 5 

4.00 

X 

10“ 5 

3.00 

X 

10“ 5 

2.00 

X 

10“ 5 

Cl 6a 

2.00 

X 

10“ 4 

1.40 

X 

10“ 4 

1.10 

X 

10“ 4 

1.00 

X 

10“ 4 

Clr 

9.09 

X 

10“ 2 

5.62 

X 

10“ 2 

3.92 

X 

10“ 2 

3.11 

X 

10“ 2 

Cl p 

-7.93 

X 

10“ 2 

-5.23 

X 

10“ 2 

-3.79 

X 

10“ 2 

-2.90 

X 

10“ 2 

Cup 

6.90 

X 

10“ 4 

-5.30 

X 

10“ 4 

-1.04 

X 

10“ 3 

-1.25 

X 

10“ 3 

Cn 6r 

-3.30 

X 

10 -4 

-1.80 

X 

10 -4 

-1.30 

X 

10 -4 

-1.00 

X 

10 -4 

C nr 

-6.90 

X 

lO -1 

-4.91 

X 

lO -1 

-4.10 

X 

lO -1 

-3.78 

X 

lO -1 

Cn p 

8.77 

X 

10“ 2 

5.42 

X 

10“ 2 

3.82 

X 

10“ 2 

3.07 

X 

10“ 2 


Table III. Aerodynamic Lift and Drag Coefficients for Design Example 



Mach 6 

Mach 10 

Mach 15 

Mach 20 

^L,a 0 

-1.329 x 10“ 3 

-1.215 x 10“ 3 

-1.251 x 10“ 3 

-1.285 x 10“ 3 

C L , a 

1.359 x 10“ 2 

1.044 x 10“ 2 

9.200 x 10“ 3 

8.700 x 10“ 3 

Cd,*o 

9.931 x 10“ 3 

6.401 x 10“ 3 

3.767 x 10“ 3 

3.325 x 10“ 3 

Co,a 

6.640 x 10“ 5 

-5.096 x 10“ 5 

-2.211 x 10“ 6 

-2.187 x 10“ 4 

C D ,a 2 

2.695 x 10“ 4 

2.405 x 10“ 4 

2.449 x 10“ 4 

2.620 x 10“ 4 
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Table IV. Trim Parameters for Four Hypersonic Cruise Conditions 


Mach 

ho, ft 

a 0 , deg 

To/W 

Vo , ft /sec 

v 2 

1 - 7%, 9 units 

6 

71000 

3.51 

0.278 

5466 

0.956 

10 

95 000 

3.83 

.217 

9 626 

.863 

15 

114 000 

2.92 

.151 

15 094 

.663 

20 

130 000 

1. 18 

.106 

22 097 

.278 


Table V. Frequency and Damping Parameters for Design Example 


Equation 

Parameter 

Control 

Parameter 

value 

Period of 
oscillation, sec 

(A26a), (A26b) 

Ca 

A 

Aileron 

Aileron 

0.70 

3.0 rad/sec 

2.1 

(A12a), (A12b) 

Cr 

LO r 

Rudder 

Rudder 

0.70 

3.0 rad/sec 

2.1 

(C7) 

<*t 

Throttle 

2.0 rad/sec 

3.1 

(5a), (5b), (5c) 

Oh 

Altitude 
outer loop 

0.1047 rad/sec 

60 

(5a), (5b), (5c) 

Ch 

Altitude 
outer loop 

0.80 




Altitude 
outer loop 

0.0698 rad/sec 

90 

(5a), (5b), (5c) 

(Ty 

Cross-range 
outer loop 

0.0785 rad/sec 

80 

(5a), (5b), (5c) 


Cross-range 
outer loop 

0.80 



Uy 

Cross-range 
outer loop 

0.698 rad/sec 

90 
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Table VI. Summary of Response Time Histories 


Figure 

Flight conditions 
(Mach number) 

Input type 

10 

6, 10, 

15, 20 

Command: 2000-ft altitude increase 

11 



Command: 20 000-ft cross-range change 

12 



Command: 20 000-ft cross-range change, polar plot of load 




factor versus bank angle 

13 



Command: combined 2000-ft altitude increase and 20 000-ft 




cross-range change 

14 

' 


Command: 10° heading change 

15 

15 

Command: 5000-ft descent, bank-to-dive versus pure descent 

16 

6, 10, 

15, 20 

Command: 5000-ft descent, bank-to-dive 

17 



Disturbance: 30-knot headwind increase 

18 



Disturbance: 30-knot crosswind 

19 



Disturbance: 10-percent density pulse 

20 

10 

Disturbance: 10-percent density pulse of various durations 
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Figure 1. Coordinate system and notation for body axes. 

Figure 2. Basic elements of hypersonic control system. 

Figure 3. Basic altitude loop proportional-integral-derivative (p-i-d) architecture. 

Figure 4. Altitude and cross-range loops interconnected through resolver. 

Figure 5. Resolver axis system and coordinate transformation. 

Figure 6. Load factor limits in resolver coordinate system. 

Figure 7. Load factor limits to prioritize altitude regulation over cross range. 

Figure 8. Maximum bank angle used in turn-rate limit. 

Figure 9. Flypersonic winged-cone configuration. 

Figure fO. Simulated perturbation responses to 2000-ft altitude command change. 

Figure fO. Concluded. 

Figure ff. Simulated perturbation responses to 20 000-ft cross-range command change. 

Figure ff. Continued. 

Figure ff. Continued. 

Figure ff. Concluded. 

Figure 12. Polar plot of load factor versus bank angle time histories for 20 000-ft cross-range command change. 
Figure 13. Simulated responses to combination of altitude and cross-range commands. 

Figure 13. Continued. 

Figure 13. Continued. 

Figure 13. Concluded. 

Figure 14. Simulated responses to a 10° heading-change command. 

Figure 14. Continued. 

Figure 14. Concluded. 
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Figure 

15. 

Figure 

15. 

Figure 

15. 

Figure 

15. 

Figure 

16. 

Figure 

16. 

Figure 

17. 

Figure 

17. 

Figure 

18. 

Figure 

19. 

Figure 

19. 

Figure 

19. 

Figure 

20. 

Figure 

20. 

Figure 

20. 


Bank-to-dive maneuver compared with pure descent at Mach 15. 

Continued. 

Continued. 

Concluded. 

Responses to bank-to-dive maneuver at four hypersonic flight conditions. 

Concluded. 

Simulated perturbation responses to 30-knot headwind. 

Concluded. 

Cross-range deviation and sideslip responses to 30-knot crosswind. 

Perturbation responses to negative 10-percent density pulse, 10 miles in length. 

Continued. 

Concluded. 

Perturbation responses to negative 10-percent density pulses of various lengths at Mach 10. 
Continued. 


Concluded. 
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